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Abstract 



We report a measurement of the differential cross section for W boson produc- 
tion as a function of its transverse momentum in proton-antiproton collisions at 
\/s = 1.8 TeV. The data were collected by the D0 experiment at the Fermilab 
Tevatron Collider during 1994-1995 and correspond to an integrated luminosity of 
85 pb" 1 . The results are in good agreement with quantum chromodynamics over 
the entire range of transverse momentum. 



Key words: PACS numbers 12.35.Qk,14.70.Fm,12.38.Qk 



Measurement of the differential cross section for W boson production provides 
an important test of our understanding of quantum chromodynamics (QCD). 
Its implications range from impact on the precision determination of the W 
boson mass to background estimates for new physics phenomena. Data from 
the production of W and Z bosons at hadron colliders also provide bounds 
on parametrizations used to describe the nonperturbative regime of QCD pro- 
cesses. 

The production of W bosons at the Fermilab Tevatron proton-antiproton col- 
lider proceeds predominantly via quark-antiquark annihilation. In the QCD 
description of the production mechanism, the W boson acquires transverse 
momentum by recoiling against additional gluons or quarks, which at first 
order originate from the processes qq' — > Wg and qg — > Wq'. When the trans- 
verse momentum ip?) and the invariant mass (M w ) of the W boson are of 
the same order, the production rate can be calculated perturbatively order by 
order in the strong coupling constant a s [1]. For p? <C M w , the calculation 
is dominated by large logarithms ~ a s ln(Mw/Pr) 2 , which are related to the 
presence of soft and collinear gluon radiation. Therefore, at sufficiently small 
Pr , fixed-order perturbation theory breaks down and the logarithms must be 
resummed [2]. The resummation can be carried out in transverse momentum 
(pr) space [3] or in impact parameter (6) space [4] via a Fourier transform. 
Differences between the two formalisms are discussed in Ref. [5]. 

Although resummation extends the perturbative calculation to lower values of 
Pr , a more fundamental barrier is encountered when p^ approaches Aqcd, the 
scale characterizing QCD processes. The strong coupling constant a s becomes 
large and the perturbative calculation is no longer reliable. The problem is 
circumvented by using a cutoff value and by introducing an additional function 
that parametrizes the nonperturbative effects [6,7]. The specific form of this 
function and the particular choices for the nonperturbative parameters have 
to be adjusted to give the best possible description of the data. 
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We report a new measurement [8] of the inclusive differential cross section 
for W boson production in the electron channel as a function of transverse 
momentum. We use 85 pb -1 of data recorded with the D0 detector during 
the 1994-1995 run of the Fermilab Tevatron pp collider. We have a ten- fold 
increase in the number of W boson candidates with respect to our previous 
measurement [9], reflecting the larger data set and an increase in electron 
rapidity coverage. An improved electron identification technique reduces the 
background for central rapidities and high pY^ by a factor of five compared 
to Ref. [9], and keeps the background contamination at a low level for large 
rapidities. Furthermore, corrections for detector resolution now enable direct 
comparison with theory. 

Electrons are detected in an electromagnetic (EM) calorimeter which has a 
fractional energy resolution of ~ 15%/ ' \J ' E[GeV) and a segmentation of Ar/ x 
A0 = 0.1 x 0.1 in pseudorapidity (rj) and azimuth (0). The D0 detector 
and the methods used to select W — > ev events are discussed in detail in 
Refs. [10] and [11] respectively. Below, we briefly describe the main selection 
requirements. 

Electron candidates are identified as isolated clusters of energy in the EM 
calorimeter that have a matching track in one of the drift chambers. In event 
reconstruction, electron identification is based on a likelihood technique [12]. 
The electron likelihood is constructed from: {%) a x 2 based on a covariance 
matrix that determines the consistency of the cluster in the calorimeter with 
the expected shape of an electron shower, (ii) the "electromagnetic energy 
fraction," defined as the ratio of the portion of the energy of the cluster found 
in the EM calorimeter to its total energy, (Hi) a measure of the consistency 
between the track position and the centroid of the cluster, and (iv) the ioniza- 
tion energy loss along the track. To a good approximation, these four variables 
are independent of each other. Electron candidates are accepted either in the 
central region, |?7 det | < 1.1, or in the forward region, 1.5 < | ?7 dct | < 2.5, where 
r/det refers to the value of 77 obtained by assuming that the particle originates 
from the geometrical center of the D0 detector. 

Neutrinos do not interact in the detector and thereby create an apparent 
momentum imbalance. For each event, the missing transverse energy ($t)i 
obtained from the vectorial sum of the transverse energy of all calorimeter 
cells, is attributed to the neutrino. 

Candidates for the W — > ev event sample are required to have an electron 
with > 25 GeV and Jfi T > 25 GeV. Additionally, events containing a 
second electron are rejected if the dielectron invariant mass M ee is close to 
that of the Z boson (75 GeV/c 2 < M ee < 105 GeV/c 2 ). A total of 50,486 
events passes this selection. 
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A major source of background stems from jets and direct photons passing our 
electron selection criteria. A multijet event can be misinterpreted as a W — > ev 
decay if one of the jets mimics an electron and there is sufficient mismeasure- 
ment of energy to produce significant $ T . The fraction of background events 
due to multijet, b quark, and direct-photon sources, also referred to as QCD 
background, is calculated by studying the electron likelihood in both a back- 
ground sample and a signal sample, as described in Ref. [13]. The total QCD 
background in the data sample is 2%; its shape is determined by repeating 
the background calculation for each p?£ bin. 

Other sources of background in the W — > ev sample are W — » tv, Z — > ee, 
and tt events. The process W — > tv — > evvv is indistinguishable from the 
signal on an event-by-event basis. To estimate this background, W — > tv 
events are generated with the same W boson production and decay model 
used in the calculation of the acceptance (see below), and the r leptons are 
forced to decay to electrons. Since the three-body decay of the r leads to 
a very soft electron p? spectrum compared to that from W — > ev events, 
the kinematic requirements keep this background to a moderate 2%. This is 
accounted for by making a correction to the acceptance for W bosons [13]. 
A Z — > ee event can be misidentified when one of the two electrons escapes 
detection or is poorly reconstructed in the detector and thereby simulates 
the presence of a neutrino. This background (0.5%) is estimated by applying 
the selection criteria to a sample of Monte Carlo Z — > ee events that were 
generated with isajet [14], processed through a GEANT-based [15] simulation 
of the D0 detector, and overlaid with events from random pp crossings that 
follow the luminosity profile of the data. The background from top quarks 
decaying to W bosons (0.1%) is estimated using herwig [16] Monte Carlo tt 
events and GEANT detector simulation. 

Trigger and selection efficiencies are determined using Z — > ee data in which 
one of the electrons satisfies the trigger and selection criteria, and the second 
electron provides an unbiased sample to measure the efficiencies. Due to the 
limited statistics of the Z — > ee data sample, we determine the shape of 
the efficiency as a function of transverse momentum using Z — > ee events 
generated with herwig, processed with a geant detector simulation, and 
overlaid with randomly selected minimum-bias pp collisions. This procedure 
models the effects of the underlying event and of jet activity on the selection 
of electrons. The efficiency for both the electron identification and the trigger 
requirements is (55.3 ± 2.2)%. 

The data are corrected for kinematic and geometric acceptance and detector 
resolution, as determined from a Monte Carlo program originally developed 
for measuring the mass of the W boson [17]. The method is described in detail 
in Ref. [18]. The program first generates W bosons with rj and p^ values cho- 
sen randomly from a double differential cross section d 2 a / dpYp drj provided as 
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input. The response of the detector and the effects of geometric and kinematic 
selection criteria are introduced at the next stage. For the present analysis, the 
input d 2 a / dpYp drj distribution is obtained using the iterative unfolding method 
described in Ref. [19]. The uncertainty due to this input distribution is eval- 
uated by using an initial distribution uniform in p^ and 77. The systematic 
smearing uncertainty is determined by varying the detector resolution param- 
eters by ±1 standard deviation from the nominal values. The total correction 
for kinematic and geometric acceptance and detector resolution for W — > eu 
events is (47.6 ± 0.3)%. 

The results for da(W — > ev)jdp^ , corrected for detector acceptance and reso- 
lution, are shown in Table 1 and plotted in Fig. 1, where the data are compared 
to the combined QCD perturbative and resummed calculation in 6-space, com- 
puted with published values of the nonperturbative parameters [6]. The error 
bars on the data points correspond to their statistical uncertainties. The frac- 
tional systematic uncertainty is shown as a band in the lower portion of the 
plot. The largest contributions to the systematic error are from uncertainties 
in the hadronic energy scale and resolution, the selection efficiency, and the 
background (in the high p^ region). An additional normalization uncertainty 
of ±4.4% from the integrated luminosity is not included in any of the plots 
nor in the table. The data are normalized to the measured W — > eu cross 
section (2310 pb [13]). The points are plotted at the values of p? where the 
predicted function equals its mean over the bin [20]. 

Figure 2 shows a comparison of the differential cross section for W boson 
production, assuming B(W — * eu) = 0.111, to the fixed-order perturbative 
calculation and to three different resummation calculations in the low p? 
region. The parametrizations of the nonperturbative region are from Arnold- 
Kauffman [5] and Ladinsky-Yuan [6] in 6-space, and Ellis- Veseli [7] in p T -space. 
The disagreement between the data and the fixed-order prediction at low val- 
ues of p^ confirms the presence of contributions from soft gluon emission, 
which are accounted for in the resummation formalisms. The fractional differ- 
ences (Data — Theory) /Theory are also shown in Fig. 2 for each of the three 
resummation predictions. Although the x 2 f° r the Ellis- Veseli and Arnold- 
Kauffman prescriptions are not as good as for Ladinsky-Yuan, the flexibility 
in parameter space and in the form of the nonperturbative function in all three 
resummed models is such that a good description of our measurement can be 
achieved [18,21]. 

Figure 3 shows the differential cross section for W boson production in the 
intermediate and high p^ regions. The calculation by Ladinsky-Yuan [6] spec- 
ifies a matching prescription which provides a smooth transition between the 
resummed and the fixed-order perturbative results to 0(a 2 ). The p^-space 
result by Ellis- Veseli [7] contains only the 0(a s ) finite part and an 0(a 2 ) Su- 
dakov form factor. Hence, there is still a residual unmatched higher-order effect 
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present in da/dp? in the large p? region, where the cancellation of the differ- 
ent parts is quite delicate. The 6-space prediction by Arnold-Kauffman [5] uses 
the matched result below p?f = 50 GeV/c and the pure conventional pertur- 
bative 0(a 2 s ) result above. We observe good agreement with the theoretical 
predictions for intermediate and high values of p^ , which probes effects of 
fixed-order QCD. 

In summary, we have used data taken with the D0 detector in pp collisions at 
y/s =1.8 TeV to measure the cross section for W — > ev events as a function 
of Pr . The combined QCD perturbative and resummed predictions are in 
agreement with the fully corrected px spectrum of W boson production in the 
kinematic range of the measurement. 
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edge support from the Department of Energy and National Science Foundation 
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Table 1 

Summary of the measurement of the px distribution of the W boson. The nominal 
Pjf is where the predicted function equals its mean value over the bin. The quantity 
da(W — > ev)/dp T 7 corresponds to the differential cross section in each bin of pjf 
for W — ► ev production. Systematic uncertainties do not include an overall 4.4% 
normalization uncertainty in integrated luminosity. 

Statistical Systematic 
T T dp^ uncertainty uncertainty 
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Fig. 1. Differential cross section for W — > ev production. The solid line is the 
theoretical prediction of Ref. [6]. Data points show only statistical uncertainties. 
The fractional systematic uncertainty, shown as the band in the lower plot, does 
not include an overall 4.4% normalization uncertainty in integrated luminosity. 
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Fig. 2. Differential cross section for W boson production compared to three resum- 
mation calculations and to the fixed-order calculation. Uncertainties on data include 
both statistical and systematical contributions (other than an overall normaliza- 
tion uncertainty in integrated luminosity). Also shown are the fractional differences 
(Data-Theory)/Theory between data and the resummed predictions. 
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Fig. 3. Differential cross section for W boson production in the intermediate and 
high pjf regions. Uncertainties on data include both statistical and systematical 
contributions (other than an overall normalization uncertainty in integrated lumi- 
nosity) . 
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